f-323 


•1- 


Ti.is  ■'  V  iJ  ort  no^c  o’.  Hjjproximits  representation  for  coiapxiting 
purposes  of  continuous  leal-valued  functions  of  several  real 
variables  by  fi  r.ite  sums  of  functions  separable  in  the  individual 
variables,  r  ) 

T 

PCLWOMI AL-LIKE  A.”  I  HOXIMATION 
Oliver  Gross 


Frequently',  in  nachine  conputations,  one  may  substitute  for  a 
function  of  relatively  complicated  analytic  structure  a  function  of 
;;impler  form  ivithout  appreciable  error,  itoreover,  if  the  function  is 
^Iven  tabularly,  a  considerable  reduction  of  tabular  input  data  (load 
on  the  memory)  may  be  achieved.  In  what  follows,  we  deal  with  xeal- 
valued  functions  of  continuous  real  variables  x,  y,  ,  Computationally, 
it  is  usually  better  to  n-or':  with  variables  assuming  finitely  many  discrete 
values.  In  tl;e  formulas  which  follow,  however,  the  discrete  analogues  are 
apparent,  integration  to  be  replaced  by  surmation,  etc,  'v'e  have  in  mind 
multivariate  approximatioi’s  by  polynomial-like  forms,  in  particular,  by 
functions  M  of  the  form 


yU»  y) 


♦  4(x)(^j(y) 
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i:i  the  livari'tte  case.  Some  of  the  folion-ing  results  are  gene’^lizable  to 
functi'.i!  9  of  more  than  two  variables. 

Of  special  interest  is  the  so-called  slide  rule  fom,  f(x)  ♦  £(y)» 
he  sl'.ciil  give  formulas  for  approjiimating  to  a  given  real-valuei  ccrtinuous 


t 


function,  defined  over  the  unit  square  0  <  ^  <  to  mlnlnize,  respectively- 
each  of  the  following  measures  of  goodness  of  fit: 


1  1 

(a)  f  y’[7(x,  y)  -  f(x)  -  g(y)J  ^dxdy,  and 

(b)  maa  |z(x,  y)  -  f(x)  -  g(y)  |  . 


A  simple  variational  technique  yields  for  (a): 
1 

-  /  z(x,  y)di'  ♦  c^  ,  and 
0 

1 

&(y)  ■  J  *(x,  y)dx  ♦  c^  ,  where 


c  and  c,  are  any  real  numbers  such  that 
0  1  *' 


®o  *  'i  ■  -  y .  • 

0  0 

Recursion  formiLlas  for  (b)  are  based  heuri^tirall^  on  the  .^^ct 
that  given  two  numbers  i  jrvl  b,  one  can  app^nKimat-:'  to  t’ -'-n  :y  a  s*  ■[I''' 
number  c;  taking  c  to  be  the  irithmetic  mean  of  a  :.nr'  h  minimi .-oa  ^  e 
maximum  absolute  error.  The  rec  jrsiou  formlar  for  (b)  arc  ac  rolJo’vc: 
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3n(y) 
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[}(■<»  y)  -  j 

[2(x,  y)  -  f^(x)]  j 


t 
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and  f  arbitrary  but  continuous,  n  •  0,  1,  ...  .  DILIPrP.TO  ♦  3T?.AU3^ 
h''.ve  shown  quite  generally  th't  the  above  process  converges  to  a  pair 
of  continuous  fxjuctions  (f,  g)  which  hao  the  require!  property.  Included 
in  their  Jiscusaion  is  a  precise  method,  discovered  independently  by  the 
.■•rlter,  for  estinating  the  value  of  the  ninimun.  In  this  connection,  the 
readei-  is  referred  to  their  paper. 

V.'e  now  return  to  the  r.ore  general  form  (1),  Lot  n  v  0  be  a  pre¬ 
scribe!  integer,  V.’e  assume  t!.'\t  z  (continuous  over  the  unit  square)  is 
not  0-'  the  form  (1)  for  this  particular  n,  V.’e  shall  give  two  sets  of 
formulas  for  approximating  to  z  by  continuous  fuiction:  of  the  ♦*orm  (1), 
.}lch  satisfy  respectively; 


(c)  The  error  vanishes  on  some  rectangulnr  grid  of  lines  x 
■p,  ...  »  "  ^1'  '•‘I*  •••  '  ^*i  ^  ^i  ^  if  i  /  J) 

y)’ -  ^  (y)j^dxdy  is  ninimizei. 


For  (c),  functions  4,  can  be  computed  successively'  an  follows, 

J 

3ho)re  a  point  (x^,  y^^)  such  that  2(3^,  y^^)  ^  0  and  two  constants  c^^,  c^ 


sue'-,  tliat  c,c. 


1  ii) 


and  set 
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-  Cj^z(x,  and 

4^(y)  "  y'*  • 

■.\«  compute  ^2  4*2  precisely  the  aaine  manner,  except  th^t  z  is  now 

replaced  by  a  new  function  22  given  by 

Z2(^*  y)  "  y)  -  ♦i(^t)4'i(y)  • 

Thus,  choose  a  point  (x^,  y2)  S'ich  that  22^*2'  ^  constants 

c«,  c^*  such  that  c^c^'  «  ~  '  r  and  set 

2  2  2  2  22(X2,  y2) 

♦2(^)  ■  ®2*2^*'  ^2^ 

4'2(y)  "  °2*®2^^?'  • 

We  continue  in  this  manner,  with 

y)  •  y)  -  (Jj^W^T^Cy) 

until  we  have  the  required  number  of  functions,  A  simple  inductive  »sr;:’ment 
shows  thit  condition  (c)  is  satisfied. 

That  the  solution  of  problem  (d),  i.e,  finding  a  least  squares 
fit  to  of  the  form  (1),  bears  a  direct  relations’nlp  to  the  theorj'  of 
integral  operators  in  Hilbert  space  was  pointed  out  by  George  Prown  at 
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RAND  a  few  years  .'ijo,  Kls  observations  are  esrontially  these; 

Let  •••  *  X.^  be  the  largest  n  eigenvalu'.s  of  the 

s^Timetric  kernel 

1 

K(x,  y)  ■/ 2(x,  t)z(y,  t)dlt 


and  •••  y  be  corresponiing  eigenfunctions  (nomalizei  so  that 

^  *  X  ’  n 

/’ (J)  .‘"(x)dx  ■  1)  and  set 
0  ^ 

1 

<l4(y)  •  f  ,  j  -  1,  2,  ...  ,  n, 

J  0 

then  the  (^*t  and  4(8  defined  in  this  wi^  fora  a  eol^jition  to  (d),  the  ^'a 

being  mutually  orthogonal  and  the  4*s  I  ing  mutually  orthogonal.  The 
com;.utatlons  can  be  carried  on  sequentially.  Let  be  the  largest  eigen¬ 
value  of  Kf  a.  corresponding  norraalizeci  eigenfiinction  and  defined  as 
above.  Having  found  (|lj^  and  4^  wo  form  the  syranetric  Kernel,  associated 
with  the  residual  function  2(x,  y)  -  •  "Ph-en,  the  largest  eigen¬ 

value  of  is  the  next  largest  eigenvalue  of  K  and  we  take  to  be  the 
corresponding  normalized  eigenfunction  wit;.  4^  defined  as  above,  V.'e  con- 
tinue  in  tlds  manner  until  we  h'lve  the  required  number  of  functions. 
Incidentally,  for  n  ■  1,  a  variational  tecluiique  yields  the  necessnry 

conditions  for  an  extreaun; 

1 

/z{x,  y)4(y)dy 

^(x)  -  ^ - 

■  2 

y  (y)^y 

o 
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<;(y) 


/z(xi  y)<>(x)<ix 
o 

— 77- - 


Ti"cse  can  ^  ujed  to  generate  an  iterative  computation  for  ^  and  (j,, 

;io\;ever|  questioiis  of  convergence,  proper  initiating  functional  to  achieve 

the  larger.t  eigenvalue,  etc.  seen  to  be  difficult.  One  fin:*!  renark  — 

the  nininun  value  of  /ylzCx,  y)  -  t^(x)t,.(y)  rdxJy  is  precisely 

00  ^  ^ 

the  sura  of  t;.e  remaining  eigenvalues  of  K,  i.e,  the  sun  of  tht  eigenvalues 
minus  the  sun  of  the  largest  n  eigenvalues,  the  former  sun  being  equal  to 
^z^dxdy. 
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